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ABSTRACT 


A study of the Southeast Asian moles of the genus Euroscaptor based on a combined approach, viz. DNA sequence 
data combined with a multivariate analysis of cranial characters, has revealed a high cryptic diversity of the group. 
An analysis of mitochondrial cytochrome b gene and five nuclear genes has revealed two deeply divergent clades: 
the western one (E. klossi + E. malayana + E. longirostris from Sichuan + Euroscaptor spp. from northern Vietnam 
and Yunnan, China), and the eastern one (E. parvidens s.l. + E. subanura). The pattern of genetic variation in the 
genus Euroscaptor discovered in the present study provides support for the existence of several cryptic lineages 
that could be treated as distinct species based on their genetic and morphological distinctness and geographical 
distribution. The moles from southern China and northern Vietnam form three distinct groups. The specimens 
from Sichuan (including the one collected from the type locality of E. longirostris) were clearly distinct from the 
northwestern Vietnam and Yunnan samples that were previously attributed to this species. We argue that the real 
distribution of E. longirostris is restricted to Sichuan, northward of Yangtze River, whereas the populations occur- 
ring southward of this isolation barrier evidently represent a new species Euroscaptor orlovi sp. nov. (northwestern 
Vietnam and Yunnan, southern China). Moreover, Red River that divides the western and eastern parts of northern 
Vietnam beyond doubt separates the population of E. orlovi from the moles occurring in northeastern Vietnam 
(Vinh Phuc and Cao Bang provinces); the latter are described here as Euroscaptor kuznetsovi sp. nov. Yet, genetic 
data are in favour of a close affinity of E. subanura with E. parvidens. A combined analysis of both genetic and mor- 
phological data has revealed a strong geographic segregation of E. parvidens samples. The populations from Dalat 
Plateau (southern Vietnam), including the moles from Loc Bao, Bi Dup and Chu Yang Sin, form a well-supported 
clade and can be considered true E. parvidens. The specimens from central Vietnam (Kon Tum and Quang Nam 
provinces) are significantly different from them, yet their monophyly has been supported by the mtDNA only. The 
moles from central Vietnam have been described here as a new subspecies Euroscaptor parvidens ngoclinhensis ssp. 
nov. All the studied samples of E. subanura have shown a low genetic and morphological variability despite their 


wide geographic range. 
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РЕЗЮМЕ 


Исследование азиатских кротов рода Euroscaptor c применением комплексного подхода, включающего ана- 
лиз ДНК и многомерный анализ краниологических признаков, выявило высокую криптическую изменчи- 
вость этой группы. Анализ митохондриального гена цитохрома b и пяти ядерных генов показал существова- 
ние двух сильно дивергировавших клад: западная клада (Е. klossi + Е. malayana + E. longirostris из Сычуани + 
Euroscaptor spp. из северного Вьетнама и Юньнани, Китай) и восточная клада (Е. parvidens s.l. + E. subanura). 
Картина генетической изменчивости рода Euroscaptor, выявленная в нашем исследовании, свидетельствует 
о существовании нескольких криптических линий, которые, на основании их генетических и морфологиче- 
ских различий и географического распространения, могут рассматриваться как отдельные виды. Кроты из 
южного Китая и северного Въетнама образуют три отдельные группы. Экземпляры из Сычуани (включая 
экземпляр из типового местонахождения Е. longirostris) заметно отличаются от кротов из северо-западного 
Вьетнама и Юньнани, которых ранее относили к этому виду. Мы считаем, что распространение E. longirostris 
ограничивается Сычуанью, к северу от р. Янцзы, тогда как популяции, обитающие к югу от этого изолиру- 
ющего барьера, относятся к новому виду Euroscaptor orlovi sp. nov. (северо-западный Вьетнам и Юньнань, 
южный Китай). Река Красная, разделяющая западные и восточные районы северного Вьетнама, отделяет 
популяции Е. orlovi от кротов, обитающих в северо-восточном Вьетнаме (провинции Виньфук и Каобанг), 
которые описаны здесь как новый вид Euroscaptor kuznetsovi sp. nov. Генетические данные свидетельствуют о 
близком родстве Е. subanura и Е. paroidens. Анализ генетических и морфологических данных выявил значи- 
тельные географические различия между исследованными выборками Е. parvidens. Популяции с Далатского 
плато (южный Вьетнам), включая кротов из Локбао, Бидупа и Чуянгсина, образуют отдельную кладу с 
хорошей поддержкой и могут рассматриваться как собственно E. parvidens. Экземпляры из центрального 
Вьетнама (провинции Контум и Куангнам) заметно отличаются от них, хотя монофилия этой группы под- 
держана только данными по митохондриальной ДНК. Кроты из центрального Въетнама описаны здесь как 
новый подвид Euroscaptor parvidens ngoclinhensis ssp. nov. Исследованные выборки Е. subanura характеризу- 
ются низкой генетической и морфологической изменчивостью, несмотря на широкое распространение. 


Ключевые слова: криптические виды, Euroscaptor, мультилокусная филогения, многомерный анализ, таксо- 
номия 
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INTRODUCTION 


Asian moles of the genus Euroscaptor Miller, 1940 
are widespread in the mountains and low hills of east- 
ern and Southeast Asia. Hutterer (2005) recognized 
six species within Euroscaptor. E. grandis, E. klossi, 
E. longirostris, E. micrura, E. parvidens and E. mizura. 
Later, morphological and genetic data (Kawada et al. 
2008) supported a species rank of the Malaysian mole 
E. malayana, which was previously included either in 
E. micrura, or in E. klossi. Recently, one more species, 
E. subanura, was described from Vietnam (Kawada 
et al. 2012). Despite the genus Euroscaptor has been 
addressed by a number of morphological and genetic 
studies (Motokawa 2004; Kawada 2005; Kawada et 
al. 2009, 2012; Zemlemerova et al. 2013; He et al. 
2014; Shinohara et al. 2015), its species composition 
and limits remain controversial. 

Molecular data by Shinohara et al. (2008) and 
Zemlemerova et al. (2013) indicated that the inclu- 
sion of the Japanese E. mizura in the genus Euroscap- 
tor had definitely made it a paraphyletic taxon. 
Furthermore, Zemlemerova et al. (2013) argued that 
E. mizura belonged to a distinct phylogenetic lineage 
of a generic rank. Further studies have corroborated 
the occurrence of a deep split between E. mizura and 
the other Asian taxa of Talpini (He et al. 2014). Un- 
doubtedly, this taxon has been awaiting its descrip- 
tion as a separate genus (see also He et al. 2014 and 
Shinohara et al. 2014). 

With the exception of E. mizura, seven species 
are now recognized in the genus Euroscaptor s.str., 
although the phylogenetic relationships among them 
remain disputable. The systematics and taxonomy of 
Euroscaptor species occurring in the Southeast Asia, 
especially in Vietnam, have not been sufficiently 
resolved yet. Vietnam bears one of the most diverse 
fauna of fossorial moles in Southeast Asia. There are 
at least five species of Talpidae including three species 
of Euroscaptor. E. longirostris, E. parvidens and E. sub- 
anura (Kawada et al. 2012; Abramov et al. 2013a). 
Previous genetic studies of Euroscaptor in Southeast 
Asia, including Vietnam (Zemlemerova et al. 2013; 
He et al. 2014; Shinohara et al. 2014), revealed a re- 
markable intraspecific variation that is indicative of 
cryptic regional species diversity. 

In the present paper, we have undertaken a multi- 
locus analysis of the phylogenetic relationships with- 
in Euroscaptor from Southeast Asia, with a special 
focus on the issues that remain controversial, such 
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as phylogenetic positions of the taxa from Vietnam. 
A multivariate analysis of cranial characters has been 
used to estimate a morphological variation among 
the Vietnamese populations. Essential changes to the 
Euroscaptor taxonomy are presented below. 


MATERIALS AND METHODS 


Taxon sampling 


The original material for the molecular phylo- 
genetic study consists of 36 specimens of three cur- 
rently recognized Euroscaptor species, 12 specimens 
of Mogera and Talpa and 16 specimens of four other 
Talpidae genera (Desmana, Galemys, Parascalops, 
Condylura). All the studied Euroscaptor specimens 
were collected in 2004—2015 from the 10 localities 
across Vietnam (Fig. 1) during the biodiversity sur- 
veys carried out by the Joint Vietnam-Russian Trop- 
ical Research and Technological Centre. Voucher 
specimens are deposited in the Zoological Institute 
of the Russian Academy of Sciences (ZIN, Saint 
Petersburg) and Zoological Museum of Lomonosov 
Moscow State University (ZMMU, Moscow). A list 
of the analysed specimens, including their collection 
localities and museum catalogue numbers, is pro- 
vided in Table 1. 

Additional 87 sequences of Euroscaptor, 10 se- 
quences of Mogera, 20 sequences of Talpa, 5 sequences 
of Scaptochirus, 18 sequences of Parascaptor and 18 
sequences of Desmana moschata, Galemys pyrenaicus, 
Parascalops breweri and Condylura cristata were 
downloaded from GenBank (Appendix 1) and used 
in the phylogenetic analyses. 


DNA extraction, PCR amplification and 
sequencing 


Pieces of muscle fixed in 96% ethanol were used 
as the sources for DNA extraction. Standard proto- 
cols for the genomic DNA extraction from ethanol- 
preserved tissues was employed, following proteinase 
K digestion, phenol-chloroform deproteinisation and 
isopropanol precipitation (Sambrook et al. 1989). 
We have sequenced the complete mitochondrial cy- 
tochrome b (cytb) gene and fragments of five nuclear 
loci: exon 11 of the breast cancer type 1 susceptibility 
protein (BRCA 1), breast cancer type 2 susceptibility 
protein (BRCA2), apolipoprotein B (ApoB), recom- 
bination activating gene 1 (КАСТ) and a-2b adren- 
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Fig. 1. Map showing the collecting localities of the studied specimens. The specimens collected during the present study are indicated 
by circles, the specimens from GenBank are indicated by squares. Locality codes see in Table 1 and Appendix 1. The type localities for 
Euroscaptor spp. are indicated by stars. 
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Table 1. List of the original material for Euroscaptor used for genetic and morphological analysis: species name, ID of specimens, collection 
and geographic location. The samples were retrieved from the following museum collections: ZMMU - Zoological Museum of Moscow 
State University, Russia; ZIN — Zoological Institute of the Russian Academy of Sciences, Saint Petersburg, Russia, BMNH - Natural 
History Museum, London, UK. 





Specimen code in 


Museum catalog 






































Species phylogenetic trees number Collecting locality and locality code (Fig. 1 ref.) 
(Figs. 2—3 ref.) where it is available 

Eupar 33 ZIN 96663 Vietnam, Kon Tum Province, Dac Glei District, 2-3 km west 
Eupar 84 ZIN 96664 f Ngoc Linh МЕ, 157057М, 107°57’E (6) 
Eupar 85 ZIN 96665 OL NEOG ININA 
Eupar AVA 14-038 ZIN 101920 
Eupar AVA 14-041 ZIN 101921 
Eupar AVA 14-074 ZIN 101922 Vietnam, Dak Lak Province, Krong Bong District, Chu Yang 
Eupar AVA 14-075 ZIN 101923 Sin National Park, 12*25/26"^N, 108°21°52°’E (7) 
Eupar AVA 14-037 ZIN 102246 

Puroscaptor paroidens Eupar AVA 14-082 ZIN 102247 
Eupar 12 ZIN 98916 
Eupar 15 ZIN 98917 Vietnam, Lam Dong Province, Lac Duong District, 5 km 
Eupar 16 ZIN 98918 NE of Long Lanh Village, Bi Dup — Nui Ba Nature Reserve, 
Eupar 17 ZIN 98919 12°10°44’’N, 108'40/44"'E (8) 





ZMMU S-173400 





Eupar 09.04.13-2 


ZMMU S-195070 











































































































Eupar 09.04.13-1 ZIN 101901 Vietnam, Lam Dong Province, Bao Lam District, Loc Bao 
Eupar 15.04.13-2 ZMMU S-195071 Forestry, 11°44.304’N, 107°42.140’E (9) 
Eupar 14.04.13-1 ZIN 101902 
ин ZIN 96318 Vietnam, Lao Cai Province, 6 km west of Sa Pa, north slope 
Eu 146 ZIN 98142 i : nq 
of Fansipan Mt., near Tram Ton forest station, 22°21°N, 
Euroscaptor sp.1 Eu 167 ZIN 9/789 103'46'E (1) 
` ZIN 101559 
BMNH 33.4.1.148 Vietnam, Lao Cai Province, Sa Pa (1) 
Eu 99335 ZIN 99335 China, Yunnan Province, Jingdong County (10) 
Eu AVA 12-288 ZIN 101531 у Я й M 
Euroscaptor sp.2 Eu AVA 12-289 ZIN 101532 гс vicinity of Tam Dao, 
ZMMU 8-161162 : 
ZN 6354 China, Sichuan Province, Chodsigou Valley (north of 
ZIN 6373 Lunganfu Town), 32°21°N, 104°20’E (11) 
BMNH 55429 
Euroscaptor longirostris BMNH 59:430 : : : 
BMNH 55.431 China, Sichuan Province, Tsao-Po 
BMNH 55.433 
BMNH 55435 
BMNH 11.2.1.24 China, Sichuan Province, Omi-San (12) 
BMNH 99.1.9 China, Sichuan Province, Yangliupa 
Eus м Vietnam, Tuyen Quang Province, Son Duong District, Hiep 
usub AVA 12-277 ZIN 101533 : 
Eusub AVA 12-278 ZIN 101534 Hoa Commune, 3 km east of Vuoc Ly Village, (northwestern 
ZIN 101535 slope of Tam Dao Mountain), 21°38°08’’N, 105°27°35’’E (3) 
Eusub AVA 13-030 ZIN 101897 
Eusub AVA 13-033 ZIN 101898 Vietnam, Ba Vi District, 50 km west of Hanoi, vicinity of Ba 
Eusub AVA 13-036 ZIN 101899 Vi Resort, 21*04/58^^N, 105°22’54’’E (5) 
Eusub AVA 13-037 ZIN 101900 
Fuüroscaptor subanura Eusub AVA 14-098 ZIN 101924 
Eusub AVA 14-103 ZIN 101925 
Eusub AVA 14-104 ZIN 101926 
Eusub AVA 14-108 ZIN 102248 
Eusub AVA 14-115 ZIN 101927 Vietnam, Phu Tho Province, Thanh Son District, Xuan Son 
Eusub AVA 14-117 National Park, 21*08^11^^N, 104°56°11”E (4) 
ZIN 102249 
ZIN 102252 
ZIN 102254 





ZIN 102255 
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ergic receptor (A2ab). The primers for amplification 
and sequencing have been retrieved from different 
publications (Teeling et al. 2000; Ohdachi et al. 2001; 
Sato et al. 2004; Shinohara et al. 2004; Lebedev et al. 
2007; Kirihara et al. 2013; Bannikova et al. 2015). 
The internal primers for A2ab sequencing have been 
designed specifically for this study: L540a 5^ -ccttgtc- 
tacctgcgcatctatctg-3' and R930 5^-gggctacaagcygag- 
gcaggaga-3'. 

The PCR protocol for all genes has been as fol- 
lows: initial denaturation at 94°С for 3 min, then 
30 cycles of 94°C for 30 s, 52—65°С (depending on 
the primer pair) for 1 min, and 72°С for 1 min, with 
a final extension of 72°C for 6 min. PCR products 
were visualised on 1.5% agarose gel and then purified 
using ammonium ethanol precipitation. Approxi- 
mately 10—30 ng of the purified PCR product have 
been used for sequencing with each primer by the 
autosequencing system ABI 3100-Avant using the 
ABIPRISMRBigDyeTM Terminator v. 3.1 (Applied 
Biosystems, Foster City, CA, USA).The sequences 
obtained in this study can be accessed via GenBank 
(Accession numbers see in Appendix 2). 


Phylogenetic analysis 


Alignment, partitioning and base composition 


All sequences have been aligned by eye using 
Bioedit version 7.0.9.0 (Hall 1999). Heterozygous 
positions (at which two peaks of approximately 
equal intensity are observed) have been coded using 
the IUB ambiguity codes. In all analyses, sequences 
were used as unphased genotypes. Phylogenetic 
reconstructions have been performed with the fol- 
lowing data sets: (1) an extended sample of taxa for 
cytb; (2) each nuclear gene separately; (2) all nuclear 
genes combined; (4) nuclear and mitochondrial cytb 
sequences combined in a species-tree estimation. 

We have determined the best partitioning strat- 
egy for nuclear concatenation using the program Par- 
tition Finder (Lanfear et al. 2012) and developed five 
a priori candidate schemes: (1) partitioning by gene; 
(2) partitioning by codon position; 3) partitioning 
by gene and codon position (three subsets per gene); 
(4) as in variant 3 but with the 1st and 2nd codon 
positions combined (2 subsets per gene) and (5) no 
partitioning. The cytb data set has always been parti- 
tioned into three codon positions. 

Pairwise genetic distances of the cytb gene se- 
quences have been estimated using the Kimura-2 
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parameter model (Kimura 1980) implemented in 
MEGA (Tamura et al. 2007). 


Phylogenetic tree reconstruction and species 
delimitation 


Phylogenetic trees have been generated by maxi- 
mum likelihood (ML), Bayesian inference (BI) and 
maximum parsimony (MP). An unweighted parsi- 
mony analysis has been performed in PAUP*4.0b10 
(Swofford 2003).The following options were in- 
voked: random addition sequence with 20 replicates, 
no limit for the number of optimal trees and TBR- 
branch swapping. Clade stability has been assessed 
based on 1000 bootstrap replicates obtained with the 
same tree search parameters. 

ML analysis has been performed in Treefinder 
(October 2008 version) (Jobb 2008). To reconstruct 
the ML tree, appropriate models of sequence evolu- 
tion selected for each partition, employing the rou- 
tine implemented in Treefinder and using BIC as the 
criterion. Separate model has been determined for 
each of the three codon positions. The following op- 
tions have been used in a tree search: parameter op- 
timisation simultaneous with tree search, optimised 
partition rates, proportional branch lengths for all 
partitions, and maximum search depth. Bootstrap 
analysis (1000 pseudoreplicates) has been performed 
with model parameters and partition rates fixed at 
the values optimal for the ML topology. 

The Bayesian analysis has been conducted using 
MrBayes 3.1.2 (Ronquist and Huelsenbeck 2003). 
Models with either two or six rate matrix parameters 
were selected for each partition based on the results 
of the model selection in Treefinder. Each analysis in- 
cluded two independent runs of four chains (one cold 
plus three heated following the default settings). The 
chain length was set at 20 million generations with 
sampling every 5000 generations. With these settings, 
the effective sample size exceeded 200 for all estimated 
parameters. Tracer 1.5 software (Rambaut and Drum- 
mond 2005) has been used to check for convergence 
and determine the necessary burn-in fraction, which 
was set to 2 million generations in all cases. 

To construct the species tree from data on six 
potentially discordant independent loci (includ- 
ing cytb), we have employed a Bayesian coalescent 
framework as implemented in *BEAST (Heled and 
Drummond 2012). The units of the analyses corre- 
spond to species or well-supported groups as inferred 
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in the Automatic Barcode Gap Discovery (ABGD) 
analysis (Puillandre et al. 2012) of the cytb alignment. 
The automatic identification of the “barcode gap" has 
been performed using the ABGD application avail- 
able at http://wwwabi.snvjussieu.fr/public/abgd/ 
abgdweb.html under the following parameters: Prior 
minimal distance 0.01, Prior maximal distance 0.1 
and relative gap width 1.0. A matrix of pairwise K2P 
distances was taken as input. 

The molecular clock assumption has been tested 
separately for each gene using likelihood ratio tests 
with likelihood values calculated in PAML 4.7 (Yang 
2007). Following the results of the molecular clock 
tests, we have used separate strict clock models for 
nuclear genes except BRCA1 and an uncorrelated 
lognormal relaxed clock model for the mtDNA and 
BRCA1. No calibration information was utilised, 
and the mean rate for A2ab was set to one. We have 
used the same partitioning scheme and models as in 
the ML analysis. Yule prior for the species tree shape 
and the piecewise constant population size model 
were assumed. Default priors were used for all other 
parameters. In total, we have conducted four runs of 
500 million generations each in BEAST version 1.8.0 
(Drummond et al. 2012). Parameter convergence has 
been assessed in Tracer 1.5. 

To validate potential species boundaries, we 
analysed the relationships between cryptic lineages 
in E. longirostris group using a multilocus species 
delimitation method, implemented in the BPP 2.1 
software (Rannala and Yang 2003; Yang and Ran- 
nala 2010). The analysis was performed based on the 
combination of one mitochondrial and five nuclear 
loci. We used Gamma (2,100) as a prior distribution 
for both 0 and «. The relative rates of the loci were 
fixed at values inferred by *BEAST. The MCMC pa- 
rameters were as follows: chain length of five million 
generations, burn-in period of ten thousand genera- 
tions, and sampling frequency set at 1/1000. 


Morphological analysis 


Eighteen cranial variables have been measured 
following Kawada et al. (2007) using a digital caliper 
with an accuracy of 0.01 mm: greatest length of the 
skull (GLS), condylobasal length (CBL), palatal 
length from the anterior tip of the 1st incisor to the 
posterior lip of the palate (PL), inner length of the 
zygomatic arch (LZA), length of the upper tooth row 
(11—M3), distance between the upper canine and 3rd 
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molar (C—M3), length of the upper molars (M 1—M3), 
rostral breadth of the canines (RB), breadth between 
the infraorbital foramina (BIOF), breadth across 
the upper 2nd molars (BAM), greatest interorbital 
breadth (IOB), braincase width (BCW), braincase 
height (BH), mandible length (MdL), mandible 
height at the coronoid process (MH), length of the 
lower toothrow (i1—m3), distance between the lower 
1st premolar and 3rd molar (p1—m3), length of the 
lower molars (m1—m3). In total, 47 skulls from the 
collections of ZIN, ZMMU and the Natural History 
Museum (BMNH, London) have been examined. 
External measurements available on museum tags, 
apparently representing the measurements obtained 
in the field by original collectors, have been used for 
comparison. 

The principal components analysis (PCA) and 
the canonical discriminant function analysis (DFA) 
have been used to evaluate a degree of cranial dif- 
ferentiation. The software program Statistica 8.0 
(StatSoft Inc., Tulsa, OK, USA) has been used for all 
analytical procedures. Males and females were com- 
bined in the analysis, because the sexual dimorphism 
was earlier reported to be insignificant (Abe 1967; 
Kawada et al. 2012). 


RESULTS 
Alignment and partitioning in the DNA analysis 


In the combined analyses of five nuclear genes, the 
final alignment consists of 4890 nucleotide positions 
including 1206 bp of BRCA 1, 951 bp of ApoB, 816 bp 
of BRCA2, 1011 bp of RAG 1 and 906 bp of A2ab. On 
the whole, the data set contains 71 specimens includ- 
ing 4 outgroups. The final alignment of cytb includes 
1140 bp for 81 specimens including 4 outgroups. 
The optimum partitioning scheme for the nuclear 
genes identified by Partition Finder under the BIC 
criterion corresponds to the variant 4 (partitioning 
by the gene and codon position but with the 1st and 
2nd positions combined). The best-fit substitution 
models employed for each of the subsets are given in 
Appendix 3. 


Phylogenetic analysis of the cytb data 
On the Fig. 2 the ML tree is shown with the val- 


ues above the branches that correspond to the Bayes- 
ian posterior probabilities in BI and the bootstrap 
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Fig. 2. The ML phylogeny of the genus Euroscaptor as inferred from the complete cytb gene sequence. Values above the branches cor- 
respond to Bayesian posterior probabilities in MrBayes (BI) and bootstrap support (1000 pseudoreplicates) in ML and MP analyses, 
correspondingly. The asterisks denote the highly supported nodes in all analyses (Bayesian posterior probabilities BPP 70.95, ML and 
MP bootstrap support 79096); filled circles mark moderately supported nodes (Bayesian posterior probabilities (BPP)>0.85, ML>70%). 
Representatives of the genera from the tribes Desmanini, Scalopini and Condylurini are used as outgroups. 
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support in ML and MP analyses. This mitochondrial 
phylogenetic tree has not supported the monophyly 
of the genus Euroscaptor, since the Japanese species 
E. mizura was recovered as the basal branch for all the 
Asian fossorial moles, while the branches of Scapto- 
chirus and one of the Parascaptor species appeared 
within the radiation of Euroscaptor. 

With the exception of E. mizura, all the Euroscap- 
tor species are subdivided into two clades. All the 
analyses (BI, ML and MP) have strictly supported a 
clade comprising Е. klossi, E. malayana, E. longirostris 
from Mouping (China, Sichuan), which is the type 
locality of this species, and an unknown form — Euro- 
scaptor sp. — from northern Vietnam and southern 
China (Yunnan). Scaptochirus moschatus shows a 
tendency to join with this clade, however the support 
of this association is low (MI, ML) or absent (MP). 
The genetic distance (K2P) between E. longirostris 
from Sichuan and Euroscaptor sp. from northern 
Vietnam and Yunnan is about 8%. Euroscaptor sp. is 
also subdivided into two subgroups: north-west one 
(Euroscaptor sp.1) and north-east one (Euroscaptor 
sp.2), with the genetic distance between them being 
6%. Within each of these subgroups, additional hap- 
logroups have been found; each corresponds to the 
specific geographic localities: Yunnan and Lao Cai 
within the first lineage, Vinh Phuc and Cao Bang — 
within the second one. 

The second clade of Euroscaptor includes E. par- 
videns s.l. and E. subanura, forming a sister group 
to that of the Parascaptor leucura branch, with 
the moderate BI and ML. The sister relationships 
E. parvidens s.l. with E. subanura shows а high support 
(K2P-distance between them - 17%). The grouping 
E. parvidens s.l. turns out to be heterogeneous: all 
geographic samples form distinct clades with K2P- 
distance between them from 6.1% to 9.5%, which in 
turn form two main subgroups: southern and central 
ones. The first one is geographically close to the type 
locality and includes the specimens from Dalat Pla- 
teau (Loc Bao, Bi Dup and Dak Lak); the second one 
originates from the Central Highlands, or Kon Tum 
Plateau (the specimens from Kon Tum, Quang Nam 
and few GenBank samples labeled as *Dak Lak", see 
Table 1 and Appendix 1). Thus, the specimens from 
Dak Lak have been allocated into two different sub- 
groups: our own sequences belong with the southern 
clade, whereas the sequences from GenBank group 
with the central clade. In the following text, we have 
named these two clades as the E. parvidens lineage 1 
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(southern Vietnam) and the E. parvidens lineage 2 
(central Vietnam: Kon Tum and Quang Nam). 


Species delimitation using ABGD and BPP 


For purposes of species delimitation, the ABGD 
method has been applied to our cytb data set (Puil- 
landre et al. 2012), which seems functional in the ab- 
sence ofreliablea priori taxonomic information. Based 
on the frequency distribution of K2P values between 
cytb haplotypes, the ABGD has retrieved 14 groups 
for primary partitions. The 14 group solution cor- 
responds to the following partition: the E. parvidens 
lineage 1, three groups within the Е. parvidens lineage 
2, two groups within Euroscaptor sp.1, Euroscaptor 
sp.2, E. longirostris, E. subanura, E. malayana, E. klossi 
and three groups within E. mizura. In general, these 
groups correspond to the main clades on the phylo- 
genetic cytb tree (Fig. 2). 

The multilocus BPP analysis provided strong 
support (posterior probability is 1.00) for the rec- 
ognition of Euroscaptor зр.1 and Euroscaptor sp.2 as 
separate lineages. 


Phylogenetic analysis of the nuclear data set 


The results of combined analysis, as well as those 
obtained on the basis of each independent gene, are 
presented in Figs 3 and 4. Neither of the alternative 
relationships obtained from individual gene trees has 
received strong support; therefore, we have analysed 
the concatenated sequence of five genes. Based on 
the nuclear concatenation, all methods have recon- 
structed identical trees with the strong support for 
most of nodes. On the Fig. 3 only the ML tree with 
the bootstrap support values for the ML and MP and 
the posterior probabilities of the Bayesian analysis on 
each node is shown. 

On the ML-tree inferred from a concatenated 
alignment of five nuclear genes (Fig. 3), the genus 
Euroscaptor (with the exception of E. müzura) is 
monophyletic in relation to Scaptochirus moschatus 
and Parascaptor, which are also monophyletic. Just 
as the case of cytb, two clades are strongly sup- 
ported within Euroscaptor s.str.: (1) E. klossi+E. 
malayana- E. longirostris from Mouping + Euroscap- 
tor spp. from northern Vietnam and southern China 
(Yunnan); (2) Е. parvidens 5.1. +E. subanura. How- 
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Fig. 3. The ML phylogeny of the genus Euroscaptor as inferred from а concatenated alignment of five nuclear genes. The designations are 
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Fig. 4. ML gene trees for Euroscaptor produced by separate analyses of the fragments of RAG1, A2ab, BRCA1, BRCA2 and ApoB. The 
asterisks denote the nodes that are highly supported in all analyses (Bayesian posterior probabilities (BPP) 70.95, ML bootstrap support 
>90%), filled circles mark moderately supported nodes (Bayesian posterior probabilities (BPP)>0.85, ML>70%). 


ever, the resolution of branches within these clades is 
not as clear as that in the mitochondrial tree. 

In the first clade, Euroscaptor sp.1 and Euroscap- 
tor sp.2 do not form a monophyletic group. The 
Vietnamese sample from Lao Cai and the specimen 
from Yunnan (Euroscaptor sp.1) are weakly grouped 


with E. malayana. The GenBank specimens of E. lon- 
girostris from Sichuan are joined with Vietnamese 
sample from Cao Bang and Vinh Phuc (Euroscaptor 
sp.2) also being weakly supported by the BI and ML 
analyses. The nuclear DNA of E. longirostris from 
Sichuan (Mouping) is presented in GenBank only by 


204 Zemlemerova et al. 


ж pEu18 Lao Саі Euroscaptor $р.1 

Eu AVA 12-288 Vinh Phuc Euroscaptor sp.2 
Eupar 14.04.13-1 Lam Dong, Loc Bao 

Eupar AVA 14-075 Dak Lak 4 
Eupar 85 Kon Tum E. parvidens 
Eupar 15 ат Dong, Bi Dup 

Eupar АУА 14-041 Рак Lak 

Eusub AVA 12-278 Tuyen Quang 

Eusub AVA 14-103 

Eusub AVA 14-108 | Phu Tho E. subanura 




















A2ab 


47 !Eusub AVA 14-104 

Eusub AVA 13-036 Ba Vi 
Mogera latouchei 
Mogera robusta 
Mogera wogura 
Talpa altaica 
Talpa europaea 
Talpa occidentalis 
Talpa caucasica 

Galemys pyrenaicus 

мм Parascalops breweri 
А Condylura cristata 


А 0.01 А 


% р НС737899 : : 
HG737900 УЕ. longirostris 
х icta LANA 228 м ян, їе " 
| Eu AVA 12-289 | Vinh Phuc> Euroscaptor sp. 
02259 | 16737901 


Eu 167 Lao Cai 
Eu 99335 Yunnan } Euroscap tor sp.1 


BRCA1 Eupar 09.04.13-1 
Eupar 09.04.13-2 |Lam Dong, 
Eupar 14.04.13-1 | Loc Bao 
Eupar 15.04.13-2 lineage | 7 ; 
2 Eupar | 1 > гапро, i Dup E. par videns 
iupar AVA 14-0 3 y 
ж ГД Eupar AVA 14-075 Jak Lak 
Eupar 12 Lam Dong, Bi Dup 
Eupar 85 Kon Tum $ lineage 2 
Eusub AVA 12-276 
Eusub AVA 12-277 
Eusub AVA 13-030 
Eusub AVA 14-103 
* Eusub AVA 14-104 » E. subanura 
Eusub AVA 13-036 
x Eusub AVA 13-037 
Eusub AVA 12-278 
Eusub AVA 14-108 
* HG737905 
HG737906 
° HG737907 
HG737908 Parascaptor leucura 
х НС737909 
х HG737910 : 
HG737911 | Scaptochirus moschatus 
DQ630251 z 
HG737898 E. mizura 
Mogera latouchei 
* 4, Mogera robusta 
Mogera wogura 
Talpa altaica 
Talpa caucasica 
Talpa europaea 
Talpa occidentalis 


{ -~—_——_o Parascalops breweri 


1.0/100 Desmana moschata 
Galemys pyrenaicus 


ondylura cristata 
0.01 
ке ————— 


Fig. 4 (continued). 
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Fig. 5. Species tree of Euroscaptor produced by the *BEAST algorithm using the Bayesian multispecies coalescent approach. Values above 
the branches correspond to Bayesian posterior probabilities. 


Table 2. Factor loadings and explained variance for the principal components PC1 and PC2 in the PCA (Fig. 6); canonical correlations 
and explained variance for the canonical axes CAN1 and CAN? in the DFA (Fig. 7). 








Characters PC1 PC2 CAN 1 CAN2 
GLS -0.907 0.276 -0.102 -0.384 
CBL -0.908 0.285 -0.098 -0.354 

PL -0.950 0.076 -0.039 -0.391 
LZA -0.387 0.793 -0.301 -0.192 
I1-M3 -0.958 -0.200 0.055 -0.422 
C-M3 -0.913 -0.285 0.081 -0.398 
M1-M3 -0.795 -0.371 0.072 -0.257 
RB -0.354 0.778 -0.203 -0.066 
BIOF -0.714 0.312 -0.046 -0.169 
BAM -0.776 0.097 0.012 -0.278 
IOB -0.482 -0.104 0.019 -0.054 
BCW -0.870 0.239 -0.090 -0.469 
BH -0.748 -0.107 -0.014 -0.332 
ML -0.930 0.222 -0.084 -0.442 
MH -0.829 0.073 -0.031 -0.296 
i1-m3 -0.915 -0.292 0.070 -0.392 
р1-т3 -0.867 -0.443 0.159 -0.455 
ш1-ш3 -0.806 -0.507 0.165 -0.345 





Explained variance (96) 64.7 13.5 55.6 28.3 
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Fig. 6. Ungrouped morphometric separation (principal components analysis) of Euroscaptor spp. 


the sequence of RAG1 (Fig. 4), which does not group 
it with the RAG1 sequences of Euroscaptor sp.1 and 
Euroscaptor sp.2. 

In the lineage of E. parvidens s.l. + E. subanura, the 
split of specimens from northern and central Vietnam 
is not as clear as in the mitochondrial tree. 

The lineage 1 of E. parvidens from southern Viet- 
nam (Lam Dong, Dak Lak) form a monophyletic 
group, however the support of lineage 2 (the samples 
from Kon Tum and Quang Nam) is absent on the 
nuclear tree. The specimens from Dak Lak are not 
united in one group: whereas our sample belongs to 
E. parvidens s.str., three specimens from GenBank 
form a separate branch. 


Species tree estimation 


The species tree inferred by using *BEAST (Fig. 
5) reproduces the topology of the five-gene nuclear 


phylogenetic tree. The deviation from the topology 
of the mitochondrial tree concerns only the low sup- 
ported branches. The relationships within Euroscap- 
tor (except for E. mizura) are again resolved in favour 
of two major clades: the highly supported sister 
relationships between E. parvidens and E. subanura 
on the one hand, and E. klossi, E. malayana, E. lon- 
girostris from Mouping and Euroscaptor spp. from 
northern Vietnam and southern China on the other 
hand. Scaptochirus moschatus forms a sister branch to 
Parascaptor leucura s.l. The basal position of E. mizura 
relative to all other moles of Southeast Asia is highly 
supported. 


Morphometric analysis 


Results of the PCA of craniometric characters in 
the revealed genetic clusters are shown in Fig. 6 and 
Table 2. As inferred from the scatterplot, E. parvidens 
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comprises three groups corresponding to the geo- 
graphical distribution of studied specimens: (1) 
Kon Tum, (2) Dak Lak, and (3) Lam Dong (Loc 
Bao and Bi Dup). The three groups diverge along 
the first principal component PC 1, in particular 
reflecting differences in the overall cranial size. The 
specimens from Kon Tum show the smallest average 
meanings of skull measurements, whereas those from 
Lam Dong - the largest one. Moles of E. longirostris 
group (sensu lato) were separated from each other 
mostly along the second principal component PC 2, 
but overlapped to some degree. The specimens from 
Sichuan regarded as E. longirostris could be diag- 
nosed by the low PC 1 and PC 2 scores in this plot. 
Within Euroscaptor sp.1 differences in PC 1 between 
samples from Lao Cai and Yunnan have been found. 
The position of Euroscaptor sp.2 was separated from 
the other two Е. longirostris group species based on a 
relatively high value of PC 2, corresponding to the 
rostral breadth and length of the zygomatic arch. 
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The morphometric differentiation in Е. subanura has 
not been revealed. 

DFAs drawing on the same variables have provid- 
ed another means of illuminating the morphometric 
distinctions between Vietnamese moles (Fig. 7, Table 
2). The discrimination between taxa of Е. longiros- 
tris group has been most strongly based upon the 
first canonical axis CAN 1, while E. subanura and 
E. parvidens lineage 2 (Kon Tum) diverged from oth- 
er taxa along the second canonical axis, in particular 
reflecting its small size. 


DISCUSSION 

Phylogenetic relationships among Euroscaptor 
Several contradicting hypotheses on phyloge- 

netic relationships among Euroscaptor species have 


been advanced. Stroganov (1948) placed longirostris 
and micrura in the separate genus Eoscalops Stro- 
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Fig. 7. Grouped morphometric separation (discriminant function analysis) drawn from the same specimens and measurements for Eu- 


roscaptor spp. as in Fig. 6. Symbols as in Fig. 6. 
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ganov, 1941, whereas other species (grandis, klossi 
and parvidens) considered in the genus Euroscaptor. 
Pavlinov (2003) distinguished two distinct species 
groups in Euroscaptor. the «longirostris»-group with 
the only species E. longirostris and the «micrura»- 
group that included all other species. Based on cra- 
nial and dental (relative length of upper molar rows) 
characters, Kawada et al. (2009) recognized two 
groups in Euroscaptor — E. grandis, E. klossi, E. ma- 
layana, E. micrura (and E. mizura) have skulls of the 
«klossi-type», whereas E. longirostris and Е. parvidens 
have skulls of the «longirostris-type». 

Contrary to these hypotheses, our data indicate 
the occurrence of two monophyletic lineages in 
Euroscaptor s.str. (without E. mizura). The eastern 
lineage includes two Vietnamese forms (parvidens 
and subanura), whereas the western lineage includes 
all other taxa (longirostris, malayana, klossi, sp.1 and 
sp.2). He et al. (2014) found out that the timing of 
divergence among these lineages could coincide 
with the Parascaptor/Scaptochirus split. A diver- 
gence between the eastern and western lineages of 
Euroscaptor-like moles seems to reflect a colonization 
history of South Asia. The widely distributed west- 
ern lineage could have rapidly dispersed throughout 
southern Asia, from China to Malaysia, in the late 
Miocene (He et al. 2014), whereas the eastern lineage 
can represent the descendants of the first radiation 
that survived in the refugees of eastern Indochina. 
Uplifting of mountains and isolation by large rivers 
could result in further speciation. However, it is nec- 
essary to study populations from India and western 
Indochina in order to infer a complete evolutionary 
history of this group, which will be a subject of a 
forthcoming study. 


Species diversity: general remarks 


The present study based on mitochondrial 
and nuclear DNA covering a wide distribution of 
Vietnamese moles has substantially extended the 
knowledge of genetic variation in Euroscaptor of this 
region. It has been found that the genus Euroscaptor 
is characterized by a high cryptic diversity. The ge- 
netic lineages which have been found for the known 
species have a clear geographical confinement. The 
essential mitochondrial genetic distances between 
these lineages (6-9%) often correspond to those 
between pairs of the sister species in other groups of 
the strictly fossorial moles such as Talpa or Mogera 
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(Tsuchiya et al. 2000; Bannikova et al. 2015). The 
phylogeographic breaks found within the Euroscap- 
tor species correlate with the definite geographical 
barriers that could block gene exchange between 
the phylogroups. It is known that moles of the genus 
Euroscaptor are distributed in the mountain areas, 
and do not occur in sympatry (Abe 1967). The only 
case is the Euroscaptor sp.2 (listed as E. longirostris) 
and E. subanura in Tam Dao, however, these species 
are really parapatric because they are separated by al- 
titudinal confinement (Kawada et al. 2012; Abramov 
et al. 2013a). The tendency to allopatric or parapat- 
ric distribution is known in other strictly fossorial 
moles of the genera Таїра or Mogera (Mitchell-Jones 
et al. 1999; Hutterer 2005). In Vietnam, different 
phylogroups of moles are associated with different 
mountain systems or separated by Red River. It is 
worth mentioning that in some cases the revealed 
intraspecific genetic forms also possess specific mor- 
phological (external and cranial) features. 

The pattern of genetic variation in the genus 
Euroscaptor revealed in the present study suggests 
the existence of several cryptic lineages that could 
be treated as distinct species based on their genetic 
and morphological distinctness and geographical 
distribution. In the following sections we shall pro- 
vide a brief discussion of the taxonomic status of the 
genetically and morphologically distinct lineages of 
Euroscaptor species occurring in Vietnam and south 
Asia and formally describe two new species and one 
new subspecies. 


Taxonomic implications 


E. longirostris lineage 


Both genetic and morphological data attest to the 
existence of cryptic species diversity within the cur- 
rently recognized E. longirostris and are in agreement 
with the data by He et al. (2014) who reported on 
the intraspecific variability between samples from 
Sichuan and northern Vietnam (Tam Dao). The 
studied specimens from Sichuan, including one from 
the type locality of E. longirostris, are clearly different 
from those from the North Vietnamese and Yunnan 
that were hitherto attributed to this species. Based 
on the orography of this region, one can assume that 
the distribution of true Е. longirostris is restricted 
to Sichuan, northward of Yangtze River. The popu- 
lations found southward of this isolation barrier 
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evidently represent a new species Euroscaptor sp.1. 
The aforementioned significant differences between 
moles from northwestern and northeastern Vietnam 
strongly corroborate their species rank. Red River 
which divides the western and eastern parts of north- 
ern Vietnam beyond doubt isolates the Euroscaptor 
sp.1 and Euroscaptor sp.2 populations. This assump- 
tion is in full agreement with the conclusion by 
Bannikova et al. (2011) who argued that Red River 
forms a distribution barrier between some Crocidura 
species. 

One of the Sichuan specimens (KIZ 0905172) 
originated from the type locality of E. longirostris 
(Baoxing, formerly Mouping), another one (KIZ 
0905290) was collected from Omei Mt., which is 
the terra typica for E. grandis Miller, 1940. The lat- 
ter was described from the single specimen (USNM 
255370) that was distinguishable by its large size 
(GLS=37 mm). The specimens from KIZ have GLS 
33.4 and 30.7 mm, i.e. within the size range of E. 
longirostris. As argued above, a sympatric occurrence 
of Euroscaptor species seems unlikely, and thus it is 
safe to consider the taxonomic status of E. grandis 
as yet obscure; it could be either an aberrant form 
of E. longirostris, or, if not, its type locality was mis- 
takenly recorded. According to the field notebook 
of D.C. Graham who collected the type specimen 
in 1930, this large mole was purchased in Kiating 
(now Leshan). Graham fixed the locality as “Mount 
Omei, Szechuan, China, at an altitude of 5000 feet” 
based on the information obtained from local hunt- 
ers. Perhaps the real collecting locality of E. grandis 
was different. More sampling in Omei Mt. region is 
needed to resolve this matter. 


E. parvidens \ineage 


A genetic differentiation of the populations of 
E. parvidens from southern and central Vietnam is 
highly pronounced, as revealed from the mitochon- 
drial cytb gene and nuclear DNA analyses. The 
populations of Dalat Plateau, including the samples 
from Loc Bao, Bi Dup and Chu Yang Sin, form a 
well-supported clade and can be attributed to true E. 
parvidens. This species was described from Di Linh 
located 40—45 km NW of our collecting site Loc Bao 
in Bao Lam District of Lam Dong, Vietnam (Fig. 1). 
The specimens from central Vietnam (Kon Tum and 
Quang Nam) are significantly different from them, 
yet, the monophyly of this group is supported by the 
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mtDNA only. This fact makes it difficult to recognize 
this clade as a separate taxon of species rank. Mean- 
while, moles from the two lineages of Е. paroidens 
differ in their morphology (viz., body and skull sizes, 
the shape of P4). The populations of Kon Tum and 
Quang Nam are rather distinct from that of Dalat 
Plateau and undoubtedly represent a new subspecies 
that is described below. It is worth mentioning that 
the high mountains of Kon Tum Plateau including 
Ngoc Linh Mt. bear a high endemic diversity of ver- 
tebrates. Ngoc Linh Mt. is known as a locality for 
such endemic species as shrew Crocidura sokolovi 
(Jenkins et al. 2007), frogs Vibrissaphora ngoclinhen- 
sis and Gracixalus lumarius (Orlov 2005; Rowley et 
al. 2014), snake Calamaria abramovi (Orlov 2009), 
and barwing Actinodura sodangorum (Eames et al. 
1999). 

The situation in Dak Lak (Chu Yang Sin, locali- 
ties 7 and 15 in Fig. 1) where both lineages are pre- 
sented remains unclear. It is hardly possible to make 
any conclusion about possible reasons for this genetic 
diversity in Chu Yang Sin area, since we know noth- 
ing about the specimens by Shinohara et al. (2014) of 
which sequences were presented in GenBank. 

The studied genetic data (see above) support a 
close affinity of the recently described E. subanura 
and Е. parvidens, as it was argued by Kawada et al. 
(2012) on the basis of morphological analysis. Nev- 
ertheless, E. subanura well differs both morphologi- 
cally and genetically from the parvidens lineage. The 
studied samples of E. subanura show low genetic 
variability, despite their wide geographic range. 
Such a low genetic diversity may be indicative of 
recent colonization of these populations across the 
most part of their distribution range during a short 
period of time. A recent colonization of E. subanura 
in Vietnam could result from a late expansion event 
northward to the areas that were already occupied 
by members of the longirostris lineage. Most of the 
high elevation areas in this region were already oc- 
cupied by /longirostris-like moles, therefore E. sub- 
anura could expand to low elevations only. In Tam 
Dao Mt., this species is distributed at the elevations 
of 200—300 m a.s.l. and another species of the longi- 
rostris group occurs at the higher elevations of 700— 
1000 m a.s.l. (Kawada et al. 2012). However, in the 
area where no sympatic species occur, E. subanura 
can occupy a higher elevation; for example, in Ba Vi 
this species occurs at the elevations of 300—1000 m 
a.s.l. (Abramov et al. 2013b). 
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Western moles 


A taxonomic diversity of the moles from western 
Indochina and India remains poorly studied. The 
Malayan mole E. malayana and Е. klossi from Thai- 
land appear to be genetically close to the longirostris- 
lineage. Shinohara et al. (2014) estimated that the 
diversification of this group (klossi, malayana and 
longirostris s.l.) took place ca. 5.3 Mya (3.8-7.9 Mya). 

The phylogenetic relationships of E. micrura and 
E. cryptura from the western Himalayas with other 
Asian moles are still unstudied. These large-sized 
moles with very short tails are known from Nepal 
and northern India. Their distribution ranges are 
isolated from the studied mole lineages by great 
Brahmaputra and Irrawaddy Rivers, so these taxa 
are likely to represent very different genetic lineages. 
Motokawa (2004) conducted a phylogenetic analysis 
of members of the family Talpidae and came to the 
conclusion that micrura and longirostris form two 
paraphyletic groups. 


SYSTEMATICS 
Genus Euroscaptor Miller, 1940 


Type species: Talpa klossi Thomas, 1929. 

Distribution: widely distributed in South Asia, 
occurring in China, Nepal, Laos, Vietnam, Myanmar, 
Thailand, and Peninsular Malaysia. 

Composition: Euroscaptor klossi (Thomas, 1929), 
Euroscaptor longirostris (Milne-Edwards, 1870), 
Euroscaptor malayana (Chasen, 1940), Euroscap- 
tor parvidens (Miller, 1940), Euroscaptor subanura 
Kawada et al., 2012, Euroscaptor orlovi sp. nov., and 
Euroscaptor kuznetsovi sp. nov. (see descriptions 
below). Probably, also includes Euroscaptor grandis 
Miller, 1940 and Euroscaptor micrura (Hodgson, 
1841). The taxonomic status of "Talpa" cryptura 
Blyth, 1843 from Darjeeling, India remains obscure. 


Euroscaptor orlovi sp. nov. 
(Figs 8B, 10B) 


Holotype. ZIN 96318, collector’s number 18, fe- 
male, body in ethanol, skull extracted, VIETNAM, 
Lao Cai Province, 6 km W of Sa Pa, north slope of 
Fansipan Mt. area, near Tram Ton forest station, 
22°21°N, 103°46’E, altitude 2000 m а..1.; collected 
10 December 2005 by A.V. Abramov. 
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Paratypes. ZIN 98142, collector's number 146, 
sex unknown, body in ethanol, skull extracted, col- 
lected 31 May 2006 by A.V. Abramov; ZIN 97789, 
collector's number 167, male, body in ethanol, skull 
extracted, collected 30 May 2008 by A.V. Abramov 
and A.V. Shchinov; ZIN 101559, collector's number 
80, female, body in ethanol, skull extracted, collected 
7 March 2012 by A.V. Shchinov; all from the same 
locality as the holotype. 

Other material. BMNH 33.4.1.148, female, skull, 
Chapa (7 Sa Pa, Lao Cai Province, Vietnam); ZIN 
99335, sex unknown, body in ethanol, skull extracted, 
Jingdong County, Yunnan Province, China. 

Diagnosis. Large-sized mole, comparable to Ё. 
longirostris. Pelage blackish brown. Tail long and 
club-shaped. Rostral part of skull elongated and nar- 
row. Posterolingual border of P4 deeply concave. An- 
terior parts of auditory bullae flatted and straddling. 
A new species distinguished by the mitochondrial 
gene cytochrome b and five nuclear genes (ВКСА17, 
BRCA2, ApoB, RAG 1 and A2ab) (see Figs 2—4). 

Description and comparisons. Head and body 
length 115—129 mm, with a long tail, 12.6—13.7?6 
of head and body length. Tail club-shaped, evenly 
covered with long overhairs. Dorsal and ventral 
pelage blackish brown. Body elongated. Skull with 
a long and narrow rostrum (see Fig. 8B). All upper 
premolars with 2 roots, in P2 roots often adnated, P2 
smaller than P1, P3 much wider them. P1 narrow, 
without metastyle. P4 with small parastyle, talon 
of a moderate size, shorter than metastyle part, the 
posterolingual border of the tooth deeply concave 
(Fig. 10B). The coronoid process with a narrow apex. 
Condyle-angular region of mandible not oblong. 

Tail long and club-shaped unlike short, bulb- 
shaped tail in E. parvidens and E. subanura. New spe- 
cies differs from E. kuznetsovi sp.nov. in the slightly 
smaller and light-built skull with a relatively narrow 
rostral part. Angular process long, band-shaped, with 
a relatively narrow base; in E. kuznetsovi sp. nov., this 
base wider and process looks triangular. The new 
species is similar in size with E. longirostris, but dif- 
fers in having short toothrows. Larger on average 
than Е. parvidens and E. subanura, with a longer tail. 
Anterior parts of auditory bullae flatted and strad- 
dling concave in E. orlovi sp. nov., E. kuznetsovi sp. 
nov., and Е. longirostris, but it inflated and converge 
in E. parvidens and E. subanura. Posterolingual border 
of P4 deeply concave in Ё. orlovi sp. nov., E. kuznetsovi 
sp. nov. and £. longirostris, but it shallowly concave in 
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Fig. 8. Dorsal, ventral and lateral views of skull, and mandible of Euroscaptor. A — E. longirostris (ZIN 6374, China, Sichuan), B — E. orlovi 
sp.nov. (ZIN 96318, holotype, Vietnam, Sa Pa), C — E. kuznetsovi sp. nov. (ZIN 101531, holotype, Vietnam, Tam Dao). Scale Баг = 10 mm. 





Fig. 9. Dorsal, ventral and lateral views of skull, and mandible of Euroscaptor. A — E. parvidens (ZIN 101901, Vietnam, Loc Bao), B — 
E. parvidens ngoclinhensis ssp.n. (ZIN 96663, holotype, Vietnam, Ngoc Linh), C — Е. subanura (ZIN 101535, Vietnam, Vuoc Ly). Scale 
Баг = 10 mm. 
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E. parvidens and E. subanura (Fig. 10). The coronoid 
process with a narrow apex, similar to E. kuznetsovi sp. 
nov. and Е. longirostris, whereas it short, with a wide 
apex in Е. parvidens and E. subanura (see Figs 8—9). 

Etymology. The new species is named in honour 
of Dr. Nikolai L. Orlov (Zoological Institute, Russian 
Academy of Sciences, Saint Petersburg, Russia) in 
recognition of his remarkable contributions to the 
study of animals of Vietnam. 

Distribution. Found in northern Vietnam (Lao 
Cai Province, Sa Pa District) and southern China 
(Yunnan Province). It may have a wider distribu- 
tion in the highlands of northern Laos and in north- 
western Vietnam, probably westward of Red River. 

Remarks. The species is reported as Euroscaptor 
sp. 1 in the results of both genetic and morphological 
analyses. 


Euroscaptor kuznetsovi sp. nov. 
(Figs 8C, 10C) 


Type locality. 

Holotype. ZIN 101531, collector's number AVA 
12-288, male, body in ethanol, skull extracted, 
VIETNAM, Vinh Phuc Province, vicinity of Tam 
Dao, 21727/06"N, 105°38°09’’E, altitude 750 m 
а.8.1, collected 24 November 2012 by A.V. Abramov. 

Paratypes. ZIN 101532, collector's number AVA 
12-289, male, body in ethanol, skull extracted, col- 
lected 25 November 2012 by A.V. Abramov; ZMMU 
5-161162, male, skin, skull, collected 6 February 
1993 by G.V. Kuznetsov; both specimens from the 
same locality as the holotype. 

Diagnosis. Large-sized mole, comparable to 
E. longirostris and E. orlovi sp. nov. Pelage blackish 
brown. Tail long and club-shaped. Rostral part of 
skull elongated and relatively wide. Posterolingual 
border of P4 deeply concave. Anterior parts of audi- 
tory bullae flatted and straddling. A new species dis- 
tinguished by the mitochondrial gene cytochrome b 
and five nuclear genes (BRCA1, BRCA2, ApoB, RAG1 
and A2ab) (see Figs 2—4). 

Description and comparisons. Head and body 
length 132—136 mm, with a long tail, 11.4—12.5% of 
head and body length. Tail club-shaped, evenly cov- 
ered with long overhairs. Dorsal and ventral pelage 
blackish-brown. Body elongated. Skull with a long 
and narrow rostrum (Fig. 8C). All upper premolars 
with 2 roots, P1 and P2 similar in height, P3 much 
wider them. P1—P2 with metastyles. P4 with small 
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Fig. 10. Upper fourth premolar in SE Asian moles. A — Е. lon- 
girostris (ZIN 6374, China, Sichuan), B — E. orlovi sp. nov. (ZIN 
96318, holotype, Vietnam, Sa Pa), C — E. kuznetsovi sp. nov. (ZIN 
101531, holotype, Vietnam, Tam Dao), D — E. parvidens (ZIN 
101901, Vietnam, Loc Bao), E – Е. parvidens ngoclinhensis ssp.nov. 
(ZIN 96663, holotype, Vietnam, Ngoc Linh), F — E. subanura (ZIN 
101535, Vietnam, Vuoc Ly). Scale bar = 10 mm. 


parastyle, talon of a moderate size, shorter than 
metastyle part, the posterolingual border of the tooth 
deeply concave (Fig. 10C). The coronoid process 
with a narrow apex. Condyle-angular region of man- 
dible not oblong. 

Tail long and club-shaped unlike short, bulb- 
shaped tail in E. parvidens and E. subanura. New 
species differs from E. longirostris in the larger skull 
size and relatively short molar rows. Larger on aver- 
age than Е. parvidens and E. subanura. It differs from 
E. longirostris and E. orlovi sp. nov. in their relatively 
wide rostral part of the skull, however it much wider 
in E. malayana and Е. klossi. The new species has P1 
and P2 similar in height, whereas in Е. orlovi sp. nov. 
and Е. longirostris P1 larger than P2. Posterolingual 
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border of P4 deeply concave in Е. kuznetsovi sp. nov., 
E. orlovi sp. nov., and E. longirostris, but it shallowly 
concave in Е. parvidens and Е. subanura (Fig. 10). 
Anterior parts of auditory bullae flatted and strad- 
dling concave in E. kuznetsovi sp. nov., E. orlovi sp. 
nov. and Е. longirostris, but it inflated and converge 
in E. parvidens and E. subanura. The coronoid process 
is narrow, differing markedly in shape from that of 
E. parvidens and E. subanura (see Fig. 9). 

Etymology. The new species is named in honour 
of Dr. German V. Kuznetsov (A.N. Severtsov Insti- 
tute of Ecology and Evolution, Russian Academy of 
Sciences, Moscow) in recognition of his many contri- 
butions to the study of mammals of Vietnam. 

Distribution. Distributed in north-eastern 
Vietnam. It is recorded from Vinh Phuc Province 
(Tam Dao) and Cao Bang Province (Nguyen Binh 
District). Found at the elevations of 750—950 m a.s.l. 
Probably, the record from Guangxi, south-eastern 
China (Hoffman and Lunde 2008) belongs to this 
species as well. 

Remarks. The species is reported as Euroscaptor 
sp. 2 in the results of both genetic and morphological 
analyses. 


Euroscaptor parvidens ngoclinhensis subsp. nov. 
(Figs 9B, 10E) 


Holotype. ZIN 96663, collector's number 33, fe- 
male, body in ethanol, skull extracted, VIETNAM, 
Kon Tum Province, Dac Glei District, 2-3 km W of 
Ngoc Linh Mt., 15*05^N, 107°57’E, altitude 1700— 
1750 m а.5.1.; collected 26 March 2006 by A.V. Ab- 
ramov. 

Paratypes. ZIN 96664, collector's number 84, 
female, body in ethanol, skull extracted; ZIN 96665, 
collector's number 85, female, body in ethanol, skull 
extracted, collected 7 April 2006 by A.V. Abramov 
from the same locality as the holotype. 

Diagnosis. Small-sized mole, comparable only to 
E. subanura, and smaller on average in its external 
and cranial measurements than the nominotypical 
E. parvidens. The fourth upper premolar with a well 
developed parastyle. A new species distinguished by 
the mitochondrial gene cytochrome b (see Fig. 2). 

Description and comparisons. Head and body 
length 120—123 mm, with a short tail, 5.8-6.1% of 
head and body length. Tail bulb-shaped, covered with 
long overhairs. Dorsal and ventral pelage dark brown. 
Body elongated with rump greatly enlarged poste- 
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riorly. Skull (Fig. 8 B) with a long and moderately 
wide rostrum. Upper and lower toothrows short. The 
upper premolars P1 and P3 similar in height and 
have 2 roots, P2 much smaller them and has one 
root. Р1-РЗ with remarkable metastyles. P4 with 
well-developed parastyle, talon wide and similar in 
size with metastyle part, the posterolingual border of 
the tooth shallowly concave (Fig. 10). The coronoid 
process short, with a wide apex. The angular process 
well-developed, with a wide base. Condyle-angular 
region of mandible oblong. 

New subspecies differs from the nominotypical 
E. parvidens in the smaller body and skull sizes and in 
the relatively short upper and lower toothrows. The 
parastyle of P4 is weakly developed in E. parvidens. 
P1 in E. parvidens parvidens is narrow, without meta- 
style, whereas in E. parvidens ngoclinhensis ssp. nov. 
the tooth is similar in width to P3 and has metasyle. 
In Е. subanura, P4 has a well-developed additional 
tubercle located anteriorly to protocone (Fig. 10F), 
which absent from both subspecies of E. parvidens. 

Etymology. The subspecies is named after the 
Ngoc Linh Mountain in Kon Tum Province of Viet- 
nam, from where it was collected for the first time. 

Distribution. Known from the Central Highlands 
of Vietnam (Kon Tum and Quang Nam provinces). 
Probably, the moles from Gia Lai Province (Abramov 
et al. 20132) belong to this subspecies. 

Remarks. The species is reported as Euroscaptor 
parvidens lineage 2 in the results of both genetic and 
morphological analyses. 
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Appendix 1. List of the GenBank accession number 
of DNA sequences used in phylogenetic analysis. For 
SE Asian Euroscaptor the geographic information 
and locality code (Fig.1) are provided. 


Euroscaptor longirostris: 

China, Sichuan Province, Baoxing (=Mouping) and Emei 
Mt. (He et al. 2014), locality 12 in Fig. 1: HG737870, 
HG737871 [cytb]; HG737914, HG737915 [RAG1]; 
HG737899, HG737900 [BRCA1]. 


Euroscaptor sp.1: 

Vietnam, Lao Cai Province, 6 km west of Sa Pa, north slope 
of Fansipan Mt., near Tram Ton forest station, 22°21°N, 
103°46’E (Zemlemerova et al. 2013), locality 1 in Fig.1: 
KC481345, KC481346, KC481347 [cytb]; KC481369 
[RAG1]; KC481357 [BRCA1]; KP717146 [BRCA2]; 
KP717258 [ApoB]; KP995409 [A2ab]; 

Vietnam, Lao Cai Province, Sapa (Shinohara et al. 2014), 
locality 1 in Fig.1: AB823111, AB823112, AB823113 
[cytb]; AB823183, AB823184, AB823185 [RAG1]; 

China, Yunnan Province, Jingdong County, locality 10 in 
Fig.1: KC481348 [cytb]; KC481370 [RAG1]; KC481358 
[BRCA1]. 
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Euroscaptor sp.2: 

Vietnam, Vinh Phuc Province, vicinity of Tam Dao, 21°27’ 
06°°N, 105°38°09’’E, altitude 750m (Meegaskum- 
bura et al. 2007; He et al. 2014), locality 2 in Fig.1: 
AB823108, AB823109, AB823110, EU122226 [cytb]; 
AB823180, AB823181, AB823182, JO433897 [RAG1]; 

Vietnam, Cao Bang Province, Nguyen Binh (Shinohara et 
al. 2014), locality 13 in Fig.1: AB823114, AB823115, 
AB823116 [cytb]; AB823186, AB823187, AB823188 
[RAG7]. 


Euroscaptor klossi: 

Thailand, Chiang Rai, Mae Sa Long, 20°09°57.4’N, 
99°37°23.2’°E, elevation 1200m (Shinohara et al. 
2014); locality 17 in Fig.1: AB823178, AB823179 
[cytb]; AB823178, AB823179 [КАСТ]. 


Euroscaptor malayana: 

Malaysia, Pahang, Cameron Highlands, 4227 М, 101°26’E, 
elevation of 1400m (Shinohara et al. 2014), locality 
17 in Fig.1: AB185155, AB185156 [cytb]; AB185155, 
AB185156 [RAG1]. 


Euroscaptor subanura: 

Vietnam, Ba Vi District, 50 km west of Hanoi, vicinity of 
Ba Vi Resort, 21704758 Ч, 10572254" E (our data), lo- 
cality 5 in Fig.1: KP995370 [cytb]; KP995375 [RAG1]; 
KP995395 [BRCA 1]; KP995388 [BRCA2]; KP995381 
[ApoB]; KP995410 [A2ab]. 


Euroscaptor parvidens: 

Vietnam, Quang Nam Province, Dong Giang District (He 
et al. 2014; Shinohara et al. 2014), locality 14 in Fig.1: 
AB823120, AB823121, AB823122 [cytb]; AB823192, 
AB823193, AB823194 [RAG]; 

Vietnam, Dak Lak Province, Chu Yang Sin National Park, 
locality 15 in Fig.1. There is no exact locality informa- 
tion in published works (Kawada et al. 2009; Shino- 
hara et al. 2014). We assume it can be Krong Bong 
District, Cu Pui Commune, based on locality data 
retrieved from Hughes (2010): AB823117, AB823118, 
AB823119 [cytb]; AB823189, AB823190, AB823191 
[RAG1]; 

Vietnam, Kon Tum Province, Dac Glei District, 2-3 km 
west of Ngoc Linh Mt., 15°05’N, 107°57”Е (Zemle- 
merova et al. 2013), locality 6 in Fig. 1: KC481342, 
KC481343, KC481344 [cytb]; KC481368 [RAG1]; 
KC481356 [BRCA7]. 

Vietnam, Lam Dong Province, Lac Duong District, 5 km 
NE Long Lanh Village, Bi Dup — Nui Ba Nature Re- 
serve, 12°10°44’°N, 108°40°44°°E (Zemlemerova et 
al. 2013), locality 8 in Fig.1: KC481338, KC481339, 
KC481340, KC481341 [cytb]; KC481365, КС481366, 
KC481367 [RAG1]; KC481354, KC481355, [BRCA1]; 
KP717147 [BRCA2]; KP717257 [ApoB]; KP995408 
[A2ab]. 
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Euroscaptor mizura: AB037604, AB076828, AB823103, 
AB823104, AB823105, DQ630413 [cytb]; AB176543, 
AB823173,  AB823174, AB823175, AB823176, 
AB823177 [RAG1]; HG737898, DQ630251 [ВКСА1]; 
DQ630168 [ApoB]. 


Mogera robusta: KC481328 [cytb]; KC481361 [RAGT]; 
KC481351 [BRCA1]; KP717144 [BRCA2]; KP717255 
[ApoB]. 


Mogera wogura: AB037646 [cytb]; HG737919 [RAGT]; 
HG737903 [BRCA1]; AB638557 [A2ab]. 


Mogera latouchei: KP717378 [cytb]; 


Talpa europaea: KF801564 [cytb]; KP717261 [RAGT]; 
KP717157 [BRCA1]; KP717098 [BRCA2]; KP717207 
[ApoB]. 


Talpa occidentalis: KP717326 [cytb]; KP717272 [RAGT]; 
KP717161 [BRCA1]; KP717106 [BRCA2]; KP717219 
[ApoB]. 


Talpa caucasica: KP717353 [cytb]; KP717298 [RAGT]; 
KP717187 [BRCA1]; KP717127 [BRCA2]; KP717240 
[ApoB]. 
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Talpa altaica: FN640579 [cytb]; KP717312 [RAGT]; 
KP717198 [BRCA1]; KP717139 [BRCA2]; KP717251 
[ApoB]. 


Scaptochirus moschatus: HG737883, AB306502 [cytb]; 
HG737927, HG737927 [RAG1]; HG737911 [BRCA1]. 


Parascaptor leucura: HG737877, HG737878, HG737879, 
HG737880, HG737881, HG737882 [cytb]; HG737921, 
HG737922, HG737923, HG737924, HG737925, 
HG737926 [КАСТ] | HG737905,  HG737906, 
HG737907, HG737908, HG737909, HG737910 
[BRCA t]. 


Desmana moschata: AB076836 [cytb]; KP717317 [RAGT]; 
KP717203 [BRCA1]; KP717148 [BRCA2]. 


Galemys pyrenaicus: AY833419 [cytb]; AY833415 [РАСТ]; 
AY121757 [BRCA1]; AY121767 [A2ab]. 


Condylura cristata: AB076810 [cytb]; KP717319 [RAG 1]; 
KP717205 [BRCA1]; KP717150 [BRCA2]; KP717260 
[ApoB]. 


Parascalops breweri: KP717379 [cytb]; KP717318 [КАСТ]; 
KP717204 [BRCA1]|; KP717149 [BRCA2]; KP717259 
[ApoB]. 
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Appendix 3. Substitution models for the mitochondrial cytb gene and five nuclear genes employed in maxi- 
mum likelihood analysis. 








cytb J2+G HKY+G GTR+G 
ApoB HKY+G K80 
BRCA1 TN+G HKY+G 
BRCA2 HKY TVM 
RAG1 Е81+ С K80+G 


A2ab K80+G K81uf+G 





